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Abstract 

In this paper, we consider an automatic-repeat-request (ARQ) retransmission protocol signaling 
over a block-fading multiple-input, multiple-output (MIMO) channel. Unlike previous work, we allow 
for multiple fading blocks within each transmission (ARQ round), and we constrain the transmitter to 
fixed rate codes constructed over complex signal constellations. In particular, we examine the general 
case of average input-power-constrained constellations as well as the practically important case of finite 
discrete constellations. This scenario is a suitable model for practical wireless communications systems 
employing orthogonal frequency division multiplexing techniques over a MIMO ARQ channel. Two 
cases of fading dynamics are considered, namely short-term static fading where channel fading gains 
change randomly for each ARQ round, and long-term static fading where channel fading gains remain 
constant over all ARQ rounds pertaining to a given message. As our main result, we prove that for 
the block-fading MIMO ARQ channel with discrete input signal constellation satisfying a short-term 
power constraint, the optimal signal-to-noise ratio (SNR) exponent is given by a modified Singleton 
bound, relating all the system parameters. To demonstrate the practical significance of the theoretical 
analysis, we present numerical results showing that practical Singleton-bound-achieving maximum 
distance separable codes achieve the optimal SNR exponent. 
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I. Introduction 

In 1957 multi-carrier transmission was first proposed by Doelz et al. [1] as a way to increase 
data rate by transmitting multiple bits streams in parallel over multiple carriers. Originally, multi- 
carrier transmission was implemented using banks of sinusoidal generators. The use of discrete 
Fourier transforms for modulation and demodulation was first suggested by Weinstein and Ebert 
in 1971 [2], significantly reducing implementation complexity, and leading to what we now 
know as orthogonal frequency division multiplexing (OFDM). A review of the development of 
multi-carrier and OFDM systems can be found in [3]. 

Almost fifty years after the invention of multi-carrier transmission [1,2], the use of OFDM 
has been adopted for broadband wireless communications systems as a means to significantly 
increase transmission rates [4]. Standards such as IEEE 802.11 (WiFi) [5,6] and IEEE 802.16 
(WiMax) [7,8] have now been extended to include OFDM techniques. Further improvements of 
data rate and reliability are promised through the use of multiple transmit and receive antennas 
[9, 10]. Multiple-input, multiple-output (MIMO) antenna systems are now being introduced into 
the IEEE 802 standards [6,8], as well as being integral parts of fourth-generation mobile 
cellular communication systems proposals [11, 12]. In addition, adaptive coding and modulation, 
combined with automatic -repeat-request (ARQ) retransmission protocols, are becoming integral 
parts of data transmission services in the Universal Mobile Telecommunications System (UMTS) 
[13], and in WiMax [8]. 

Practical wireless communication systems will therefore soon feature MIMO OFDM modu- 
lation with overlaying ARQ protocols. It is thus important to obtain a thorough understanding 
of the fundamental characteristics of such systems. In this paper, we model a practical point- 
to-point MIMO OFDM ARQ wireless communication system as a system transmitting signals 
from a complex signal constellation over a block-fading MIMO ARQ channel. In the following 
subsections, we first review prior art and technical concepts relevant to our work. We then for- 
mulate our problem and summarize contributions, before outlining the organization and defining 
notation of the paper. 

A. Prior Art 

1) Fundamental Tradeoff: The work of Teletar [14], and Foschini and Gans [15], has inspired 
a flurry of research activities in MIMO antenna systems for wireless communications. Previously, 
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multiple-antenna systems were primarily used for providing receiver diversity, thus combatting 
random amplitude fluctuations due to fading [16]. In contrast, the prevailing thesis for MIMO 
systems is that fading can increase channel capacity by providing a set of well-behaved parallel 
channels [14, 15]. In fact, in the high signal-to-noise (SNR) regime it has been shown that the 
capacity of a channel with N t transmit antennas, N r receive antennas, and independent, identical 
distributed (i.i.d.) complex Gaussian channel gains between each antenna pair is given by 

C(SNR) = nun{/V t ,/V r .}logSNR + 0(l), 

suggesting that capacity increases linearly with the minimum number of transmit and receive 
antennas. Therefore, the use of multiple-antenna systems can improve both reliability and data 
rate, when transmitting over a quasi-static MEMO channel where channel gains are i.i.d. complex 
Gaussian and fixed during the transmission. 

Zheng and Tse described the fundamental tradeoff between diversity gain and multiplexing 
gair3 for quasi-static MIMO channels in the high SNR regime in [17], assuming Gaussian dis- 
tributed input signals. The fundamental tradeoff developed in [17] has since become a benchmark 
for the performance evaluation of space-time coding schemes, and the corresponding framework 
has become a preferred approach for characterizing classes of MIMO channels. For example, in 
[18] the fundamental diversity-multiplexing-delay tradeoff is characterized for the MIMO ARQ 
channel, and the fundamental diversity-multiplexing tradeoff for MIMO channels with resolution- 
constrained feedback is determined in [19], both under the assumption that Gaussian distributed 
input signals are used. 

2) OFDM and the Block-Fading Channel: The block-fading channel model was introduced in 
[20], with the purpose of modelling slowly varying fading channels where the duration of a block- 
fading period is determined by the channel coherence time. Within a block-fading period, the 

'The diversity gain (or SNR exponent) is defined as 

SNR^oo log SNR 

where P e (SNR) denotes the probability that the transmitted message is decoded incorrectly. The multiplexing gain is defined 
as 

A fl(SNR) 
rm - SN R^oc log SNR' 

where i?(SNR) is the code rate. The multiplexing gain essentially quantifies how close the code rate is to the capacity of a 
single-input single-output link at high SNR [17]. 
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channel fading gain remains constant, while between periods the channel gains change randomly 
according to a fading distribution. In this setting, transmission typically extends over multiple 
block-fading periods. A thorough treatment of fading channels is found in [21]. 

The block-fading channel model is a reasonable model for OFDM transmission over frequency- 
selective wireless channels, as an OFDM system is typically designed such that each sub-carrier 
experiences flat fading. Despite its simplicity, the model captures important aspects of OFDM 
modulation over frequency- selective fading channels and proves useful for developing coding 
design criteria. 

The definition of multiplexing gain, fundamental in the formulation presented in [17, 18], relies 
on coding schemes with transmission rates that increase linearly with the logarithm of the SNR. 
Non-zero multiplexing gains can only be achieved with continuous input constellations or discrete 
constellations with cardinalities scaling with the SNR. From a practical perspective, it is desirable 
to operate at a fixed code rate and deal with small alphabet sizes. We are therefore interested 
in the performance of such practical schemes, which effectively operates at zero multiplexing 
gain. Under this scenario, the general diversity-multiplexing tradeoff can only provide a coarse 
characterization of the rate-diversity tradeoff. The rate-diversity tradeoff for fixed-rate space-time 
codes constructed over discrete signal constellations, and transmitted across a quasi-static MIMO 
channel, was presented in [22]. 

Union-bound arguments [23] and error exponent calculations [24] were used to show that the 
diversity gain of a block-fading channel with an arbitrary, but fixed number of fading blocks, 
fixed code rate, and a discrete input signal constellation, is described by a modified version of 
the Singleton bound [25]. The same problem is considered in [26], where outage probability 
arguments are used to formally prove that the optimal rate-diversity tradeoff is indeed the 
modified form of the Singleton bound presented in [24,25], which is achieved using maximum 
distance separable (MDS) codes. 

The block-fading ARQ channel model has recently been considered in [27, 28] for discrete 
input signal constellations. In [27] the Singleton bound is presented as an upper bound to the 
SNR exponent, while the optimality of the Singleton bound is formally proven for the ARQ case 
in [28]. In [28] it is also demonstrated that asymptotically optimal throughput can be achieved 
by MDS codes. 
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B. Problem Formulation and Contributions 

In this paper, we consider an ARQ system signaling over a block-fading MIMO channel with 
L maximum number of allowable ARQ rounds and B fading blocks per ARQ round. In contrast 
to the work in [17, 18], we allow for multiple fading blocks within each transmission (ARQ 
round), and we constrain the transmitter to fixed rate codes constructed over complex signal 
constellations. In particular, we examine the general case of average input-power-constrained 
constellations as well as the practically important case of discrete constellations of finite car- 
dinality. The receiver is able to generate a finite number of one-bit repeat-requests, subject to 
a latency constraint, whenever an error is detected in the decoded message. A maximum of L 
transmissions pertaining to each information message is allowed. 

As in [18], we consider two cases of fading statistics; for the short-term static fading case, 
the channel fading gains change randomly for each ARQ round, while for the long-term static 
fading case, the channel fading gains remain constant over all ARQ rounds pertaining to a given 
message, but change randomly for each message and corresponding suite of ARQ rounds. This 
scenario is a suitable model for practical wireless communications systems employing OFDM 
modulation over a MIMO ARQ channel. 

The main focus of our work is to derive the optimal tradeoff between throughput, diversity 
gain and delay of ARQ schemes signaling over block-fading MIMO channels. In particular, we 
show that the tradeoff highlights the roles of the complex-plane signal constellation, the rate of 
the first ARQ round Ri, the maximum number of ARQ rounds L, and the number of fading 
blocks per ARQ round B. 

As a first result, we prove that for the block-fading MIMO ARQ channel with the input 
constellation satisfying a short-term power constraint, the optimal SNR exponent is given by 
N t N r LB for short-term static fading and N t N r B for long-term static fading, which is achieved 
by Gaussian codes of any positive rate. This is, however, not the case with discrete signal 
constellations. In order to attain full diversity the signal constellations must feature certain 
properties. In general, due to the discrete nature of these signal sets, a tradeoff between rate, 
diversity and delay arises. 

As our main result, we prove that for the block-fading MIMO ARQ channel with discrete 
input signal constellation of cardinality 2 QNt satisfying a short-term power constraint, the optimal 



February 1, 2008 



DRAFT 



6 



SNR exponent is given by a modified Singleton bound, relating all the system parameters. Note, 
however, that modulating across all fading blocks increases the dimensionality of the decoding 
problem by a factor of B [29]. For further flexibility in terms of decoding complexity, we 
consider the case where modulation is performed over a number 1 < M < B of fading blocks, 
such that B = MD. The resulting optimal SNR exponent is then expressed as 



The expression in (QQ) implies that as the target rate Ri increases, the achievable optimal diversity 
order d* D (Ri) decreases in steps. Our main result generalizes the result of [22] for the quasi-static 
MIMO channel to the ARQ block-fading case with encoding across M fading blocks. 

Directly following from the results, we demonstrate that while the optimal SNR exponent of 
the system is an increasing function of the maximum number of allowed ARQ rounds L, the 
throughput of the system becomes independent of L for sufficiently high SNR, and is determined 
by the rate of the first ARQ round. We therefore denote our main result as the optimal throughput- 
diversity-delay tradeoff. This result provides strong incentive to use ARQ as a way to increase 
reliability without suffering code rate penalties. 

To demonstrate the practical coding aspects of our results, some examples are presented with 
corresponding error rate and throughput performances. The diversity tradeoff function can be 
viewed as a modified version of the Singleton bound [25], which naturally leads us to investigate 
the role of Singleton-bound-achieving MDS codes. Our examples illustrate that the optimal SNR 
exponent can be achieved with practical MDS coding schemes. 

C. Organization and Notation 

The paper is organized as follows. In Section |n] we define the system model, and in Section 
ITTll we review relevant ARQ performance measures, namely, error probability, throughput and 
average latency. In Section [IV] we review the concepts of information accumulation and out- 
age probability, while the main theorems of the paper, detailing the throughput-diversity-delay 
tradeoff, are presented in Section [V] A thorough discussion is included in Section [V] where 
the results are interpreted and related to existing results in the literature. To demonstrate the 
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practical relevance of the results, numerical examples are included in Section [VH showing that 
MDS codes achieve the tradeoff. Concluding remarks are summarized in Section IVIIL while the 
details of the proofs have been collected in the appendices. 

The following notation is used in the paper. Sets are denoted by calligraphic fonts with 
the complement denoted by superscript c. The exponential equality f(z) = z d indicates that 
lirn^oo = d. The exponential inequality <,> are similarly defined, y and -< denote 

component-wise inequality of > and <, respectively. I denotes the identity matrix, vector/matrix 
transpose is denoted by ' (e.g. v') and || ■ \\p is the Frobenius norm. !{•} is the indicator function, 
and \x] ( [x\ ) denotes the smallest (largest) integer greater (smaller) than x. 

II. System Model 

In this section we describe the block-fading MIMO ARQ channel model and coded modulation 
schemes under consideration. 

A. Channel Model 

Consider a block-fading MIMO ARQ system with N t transmit antennas and N r receive 
antennas. We investigate the use of a simple stop-and-wait ARQ protocol where the maximum 
number of ARQ rounds is denoted by L. Each ARQ round consists of B independent block- 
fading periods, each of length T (coherence time/bandwidth) in channel uses. Hence each ARQ 
round spans BT channel uses. Figure CD shows the overall system model. We write the received 
signal at the 6th block and £th ARQ round as 



where e C NtXT ,Y ttb ,W e<b G C N ' xT and H ib e C^*^ denote the transmitted signal 
matrix, received signal matrix, the noise matrix and the channel fading gain matrix, respectively. 
We define x.g b,t ^ Nt as the vectors containing the transmitted symbols of each antenna at 
ARQ round i, block b and time t, which are such that = [x^6 1, . . . , x^y]. 

Both the elements of the channel fading gain matrix H4& and the elements of the noise 
matrix are assumed i.i.d. zero mean complex circularly symmetric complex Gaussian with 
variance a 2 = 0.5 per dimension. We assume perfect receiver-side channel state information 
(CSI), namely, the channel coefficients are assumed to be perfectly known to the receiver. 




(2) 
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We obtain the long-term static model of [18] by letting H tjb = H^ )6 for all £ ^ £' in ©, 
namely, all ARQ rounds undergo the same MIMO block-fading channel. This models well 
a slowly varying MIMO OFDM ARQ system with B subcarriers or B groups of correlated 
subcarriers. On the other hand, when the matrices H^ 6 are i.i.d. from block to block and from 
ARQ round to ARQ round, © corresponds to the short-term static model of [18]. In order to 
keep the presentation general, and since © encompasses both models, we will index the channel 
matrices according to ARQ round and block as in the short-term static model. We will outline 
the changes for the long-term static model whenever necessary. 

Therefore, the channel model corresponding to ARQ round I becomes 

+W,. (3) 

where 

= \y' i)1: • • • , Y^ jB ] G C i 

= [X^ l3 . . . , X^ B ] G C BNt/ 

H, = diag(H,, 1 ,...,H,, B )GC^ x ^. 

One channel use of the equivalent model © corresponds to BT channel uses of the real channel 
©. In a similar way to the previous model, we define the vectors x^t G C BNt for t = 1, . . . , T 
as 

X,= [x £il ,...,x £ , T ]GC BiV ' xT 

The receiver attempts to decode following the reception of an ARQ round. If the received 
codeword can be decoded, the receiver sends back a one-bit acknowledgement signal to the 
transmitter via a zero-delay and error-free feedback link. The transmission of the current code- 
word ends immediately following the acknowledgment signal and the transmission of the next 
message in the queue starts. If an error is detected in the received codeword before the Lth 
ARQ round, then the receiver requests another ARQ round by sending back a one-bit negative 
acknowledgment along the perfect feedback path. However, a decision must be made at the end 
of the Lth ARQ round regardless of whether errors are detected. 

In general, the optimal ARQ decoder makes use of all available coded blocks and correspond- 
ing channel state information up to the current ARQ round in the decoding process. This leads 

February 1, 2008 DRAFT 



9 



to the concept of information accumulation, where individual ARQ rounds are combined, along 
with any other side information. We hence introduce the ARQ channel model up to the tth ARQ 
round, completely analagous to ©, but allowing for a more concise notation. In particular, we 
have that 

Y e = ^Zh^ + W,, (4) 

where 

Y f = [Y;,...,Y;fec fMrxr , 

x^[xi,... 1 xa'GC Mtxr , 

w f = [w' 1 ,..,w;fGC™ rXT , 

H, = diag(H 1 ,...,H £ )GC^' x ^. 

That is, Yi, and are simply collections of the received, code and noise matrices, re- 
spectively, available at the end of the £th ARQ round, concatenated into block column matrices. 
The new channel matrix G C eBNrX£BNt is a block diagonal matrix with the diagonal blocks 
composed of the respective channel state during each block-fading period up to ARQ round t. 
In the case of long-term static model, = diag (H, . . . , H). Note that a channel use of the 

l times 

equivalent model (HJ) corresponds to IBT channel uses of the real channel ([2]). 
B. Encoding 

In this section we discuss the specific construction of the space-time ARQ codewords. The 
information message m to be transmitted is passed through a space-time coded modulation 
encoder with codebook C C C LBNtxT and code rate R , where i? — 7 1 and 

Ri = ^ lo S2 \c\ 

is the code rate of the first ARQ round. Therefore, \C\ = 2 RoLBT and m G A4, where M. = 
{1,2,..., 2 RoLBT } is the set of possible information messages. We denote the codeword corre- 
sponding to information message m by X(m). The rate R codeword can be partitioned into a 
sequence of LB space-time coded matrices, denoted X^j G C NtxT . According to the previously 
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described model, we have that 

X(m) = [X' 1 (m),... 1 X' L (m)}' 

= [X' y ( ffl ), . . . .xym), • • • ,X' Ljl (m), . . . ,X' i)B (m)]' G C L ^ xT 

We consider a s/zo/t term average power constraint, namely, the transmitted codewords are 
normalized in energy such that VX G C, -|-4fE[||X|||.] = AT t . Therefore, together with the 
model assumptions in the previous section, p in (O, © and © represents the average SNR per 
receive antenna. 

In this paper we analyze space-time coded modulation schemes constructed over discrete 
signal sets. In particular, we consider that C is obtained as the concatenation of a classical 
coded modulation scheme Cq C Q LB ™t constructed over a complex-plane signal set Q = 
■ • • j Q\Q\} C C [30] with a unit rate linear dispersion space-time modulator [31]. Let c Q G Cq 
denote a codeword of Cq of length LBTN t and Q = log 2 \Q\ the number of bits conveyed in 
one symbol of Q, namely, |<2| = 2 Q . Since the linear dispersion space-time modulator has unit 
rate we have that < i?i < N t QL. 

To allow for a general case, we consider that the linear dispersion space-time modulator 
spreads the symbols of c Q over the N t transmit antennas and the B fading blocks. In particular, 
we consider that the codewords c Q of Cq, of length LBTN t are partitioned into L vectors 
of length BTN t each, denoted by c Qti G Q BTNt such that c Q = [c' Q1 , . . . , c'q L ] . For every 
£ = 1, . . . , L, the vectors eg >e are multiplied by the unit rate generator matrix of the linear 
dispersion space-time modulator R G C BTNtxBTNt to form 

x^ = Rc Qi ^ (5) 

where = vec(X£) G C BNtT is the vector representation of the portion of codeword of C G Cq 
transmitted at ARQ round £. Without any loss in generality we consider that R is a rotation 
matrix [32-35], i.e., R is unitary [36]. Note that introduction of the linear dispersion space-time 
modulator rotation matrix R increases the decoding complexity compared to the unrotated case 
where R = I. This is due to the fact that now the components of depend on each other, since 
R induces a change of the reference axis for detection [32-35]. This implies that the detection 
problem is of dimension BTN t . 
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To allow for further flexibility, we consider the case where the linear dispersion space-time 
modulator spreads the symbols of eg G Cq over the N t transmit antennas and a number 1 < 
M < B of fading blocks, such that 

M 

is an integer representing the number of rotations used in an ARQ round. In this case, we have 
that the rotation matrix R becomes block-diagonal, namely 

R = diag( R M ,Rm ) (6) 

D times 

where R M G C MNtTxMNtT is the rotation matrix of dimension MN t T x MN t T. According to 
© we can define G C MNtT , such that 5q = [x^i, . . . , x^,d]'. We define the multidimensional 
constellation X M as 

X M t {x G C MN * T : Vc G Q MNtT , x = R M c} (7) 

Due to the block-diagonal structure of R, the detection problem reduces to D detection problems 
over X M each of dimension MTN t . This formulation encompasses many cases of interest, as for 
example the unrotated case, for which R = I, the general threaded algebraic space-time (TAST) 
modulation structure for MIMO block-fading channels [29], or perfect space-time modulation 
[37]. As we shall see in Section |V] the parameter M plays a key role in the reliability of 
the overall system. Intuitively, the larger M, the larger the space-time symbol spreading, and 
hence, the larger the diversity [29]. On the other hand, using large M implies larger decoding 
complexity, as the detection problem is exponential in M. Using the previous discussion, we 
introduce the following equivalent channel matrix 

= diag(H^(d_i)jvf + i, . . . , THe,(d-i)M+i, ■ ■ ■ > Hg )dM , • • • , H^m ) £ c M7VrTxM7V * T (8) 
Vv v ' s v '> 

T times T times 

for d = 1, . . . ,D. These matrices correspond to the channels seen by rotation d within ARQ 
round I. The equivalent channel defined by ([8]) induces the following channel model 

9e.d = J -rrH-e,d%,d + v?e,d (9) 
V 

where x^ G C MNtT , y^, W£ t d G C MNrT are the corresponding input, output and noise vectors. 
This model describes the relationship between the output of one of the D rotations of the 
linear dispersion space-time modulator and the output of the channel. One use of channel © 
corresponds to MT uses of the real model ©• 
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C. Decoding 

We will make use of the ARQ decoder proposed in [18], which behaves as a typical set 
decoder for the first L — 1 ARQ round and finally performs ML decoding at the last ARQ 
round. The decoding function at ARQ round £, for £ = 1, . . . , L — 1, denoted ipt(Yi, H^), gives 
the following output 

m if X(m) is the unique codeword in C jointly typical with Y £ given 



otherwise, 



(10) 



which implies that message index ip^iY^ H^) = rh G M. whenever the received matrix can be 
decoded and ^(Y^H^) = whenever errors are detected. 



III. ARQ Performance Metrics 

In this section we introduce a few performance metrics relevant to ARQ systems, namely, the 
error probability, average latency and throughput. For ease of notation, we define three relevant 
decoder events as follows. Let, 



V t ± (Vi(Yi, Hi) = 0, ... , ^(Y £ , H e ) = 
denote the event of error detection up to and including ARQ round I, let 

denote the event of decoding a valid message at ARQ round t, and let 

(J ^(Y*H*)=m 

denote the event of a decoding error at ARQ round t, given that message m was transmitted. 
Based on the events defined above, the probability of error P e (p) is given by 



Pe(p) = E 



L-1 



Pr (A u S 1 ) + y2 Pr (V t -u At, £ t ) + Pr £ L ) 



=2 



(ID 



undetected errors ML decoding errors 

where the expectation is with respect to the joint distribution of the fading gain matrix and 
received signal matrix. From the error expression in (fTTj) it is clear that the ARQ decoder suffers 
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from undetected errors and ML decoding errors. Undetected errors occur during ARQ rounds 
£ = 1 . . . L— 1 and reflects the inability of the decoder to identify erroneous frames. ML decoding 
errors occur at the last ARQ round and reflects the inability of the decoder to resolve atypical 
channel and noise realizations. We shall see later that the probability of undetected errors can 
be made arbitrarily small using appropriate codebooks, leaving ML decoding errors to dominate 
the error probability. In terms of error probability, the effectiveness of an ideal ARQ decoder is 
therefore almost exclusively limited by the error probability at the last ARQ round. 

The expected latency k of the system is determined by the probability of error detection, and 
it is given by 

L-l 

K = l + ^Vr{V e ), (12) 

t=i 

where k is expressed in terms of number of ARQ rounds. The corresponding transmit throughput 
of the system in terms of the average effective code rate is simply obtained by 

v(RuL) = _ — , (13) 

1 + E|=i Pr (A) 

where ry(i?i, L) is expressed in bits per channel usa^. 

IV. Information Accumulation and Outage Probability 

In this section, we expand on the idea of mutual information accumulation in ARQ systems 
as well as introduce the commonly used concept of information outage. 

The instantaneous input-output mutual information of the channel © up to ARQ round £, for 
the channel realization H<? = G e can be written as 



l{p\G e ) = ~I(X< ; Y, | H £ = G t ). (14) 
1 1 

= -J2l(p\G £ ) (15) 



k=\ 

where J(p|G^) is the instantaneous input-output mutual information corresponding to ARQ 
round I. Following (fl"5l) we will refer to J(p|G^) as the accumulated mutual information up to 
ARQ round £. The accumulated mutual information /(p|G^) measures the normalized mutual 

2 Note that our definition of transmit throughput here is purely a measure of the average code rate at the sender's side, as it 
does not take into account whether or not messages are correctly decoded at the receiver's side. 
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information between the accumulated received matrix and the coded blocks X^, given the 
instantaneous channel state matrix Gi. Since Gi is a random matrix, l(p\GA is a non-negative 
random variable. Further, from (fl"5l) it is clear that the accumulated mutual information is an 
increasing function of the ARQ round index £, for a given realization of G^. 

Following [38, Lemma 1], we get that for \Ai\ = 2 RlBT , there exists a codebook C such 
that the conditional probability of error P e (p\Gi) < e for any e > whenever the accumulated 
instantaneous mutual information satisfies l(p\Gi) > R± for any £ = 1, . . . ,L, provided that 
the block length £BT is sufficiently large. We hence define information outage as the event that 
occurs when the accumulated mutual information is below Ri, namely 

O t 4 {g, G C eBTN ^ eBTNt : l(p\G e )< Ri} ■ (16) 

For any finite B and L, the channel defined in © is not information stable and the channel 
capacity in the strict Shannon sense is zero [39], since the probability of the outage event is 
nonzero. The corresponding outage probability is defined as [20,21] 

P mt (p,£,R 1 )=Pr(O i ) (17) 
= Pr(/(p|G<) <R,y (18) 

The accumulated mutual information l(p\Gi), and hence the corresponding outage probability, 
depends on the SNR p and the input distribution -Px(X) with the constraint that -^=E[||X||^] = 
N t . When no other constraints are imposed on the input distribution, the input distribution that 
maximizes l(p\Ge) and therefore minimizes P out (p,£, R±) is the Gaussian distribution, namely, 
the entries of X are i.i.d. complex circularly symmetric random variables with zero mean and 
unit variance. This leads to, 

I(p\G e ) = -^log 2 det (l + A-GcGl) (19) 

= \ S log2 det C 1 + F G ^ G U • (20) 

b=i \ t / 

In practice, Gaussian codebooks are not feasible, and we will resort to discrete signal constella- 
tions. In this work, we are mostly interested in studying the role of the discrete nature of practical 
constellations, and the impact this further system constraint has on the outage probability. In 
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particular, we can write the mutual information for the scheme described in Section III-BI as, 

1 D 

I(p\Ge) = -Y, I (p\Gi4) (2D 

where 



D 

d=l 




e 



/^G £ , d (x-x')+w|| 2 + ||w|| 2 



(22) 



(23) 



is the input-output mutual information corresponding to the realization He,d = G^d given in ([8]) 
of the channel described in ©, assuming a uniform distribution over the MN t T multidimensional 
constellation Xu defined in ©. Since < I(p\G^d) < QN t , it is not difficult to show that (|2T|) 
can be bounded as follows 

1 D i 1 M ( \ 1 

/(p|G,) < - < QN t , — Y log 2 det I + J^G Ud „ 1)M+m G{ [d _ l)M+m . (24) 

d=l I m=l V i / J 

This relationship will prove useful in proving our main results. 

V. Throughput-Diversity-Delay Tradeoff 

In this section, we derive the optimal tradeoff between throughput, diversity gain and delay 
of ARQ schemes signaling over MEMO block-fading channels. In particular, we show that the 
tradeoff highlights the roles of the complex-plane signal constellation through Q, the rate of the 
first ARQ round R\, the maximum number of ARQ rounds L and the number of fading blocks 
per ARQ round B. As we shall see, for large SNR, the tradeoff expression highlights the role 
of the asymptotic throughput through R\. Furthermore, the optimal tradeoff expression includes 
the effect of the space-time spreading dimension of the linear dispersion modulator, providing 
also a reference of decoding complexity. 

We now present the main results of this paper concerning the optimal SNR exponent of ARQ 
systems. 

Theorem 1: Consider the channel model © with input constellation satisfying the short term 
average power constraint -^-|^E[||X||^] < N t . The optimal SNR exponent (i*(i?i) is given by 

N t N r LB for short-term static fading 
d*(Ri) = { (25) 

N t N r B for long-term static fading 
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Further, this is achieved by Gaussian random codes of rate Ri > 0, provided that the block 
length is sufficiently long. 

Proof: Theorem 1 follows immediately as a corollary of [18, Theorem 2] after taking into 



Theorem 1 states that Gaussian codes achieve maximal diversity gain for any positive rate. As 
we show in the following, this is not the case with discrete signal constellations Xu- In particular, 
full diversity is achievable by discrete signal sets provided the rates satisfy < R\ < QN t L. 
However, in order to attain full diversity we must restrict the signal constellations to certain 
properties. In general, due to the discrete nature of these signal sets, a tradeoff between rate, 
diversity and delay arises. This relationship is expressed in the next theorem. 

Theorem 2: Consider the channel model © satisfying the short term average power constraint 
2^E[||X|||,] < N t , with discrete input signal constellations of cardinality 2 QNt . The optimal 
SNR exponent is given by 



over the full range of < Ri < QN t L where (|26l) is continuous. 

Proof (Sketch): A sketch of the proof is provided here, with the technical details left to 
Appendix. We first prove the converse and show that the diversity gain d* D (Ri) is upper- 
bounded by (|26l) . We can use Fano's inequality to show that the outage probability P ou t(p, Ri) 
lower-bounds the error probability P e (p) for a sufficiently large block length. Then we bound 
the maximum SNR exponent by considering the diversity gain of the outage probability. For 
large SNR, the instantaneous mutual information is either zero or QN t bits per channel use, 
corresponding to when the channel is in deep fade and when the channel is not in deep fade, 
respectively [26]. Achievability is proved by bounding the error probability of the typical set 
decoder [18] for ARQ rounds I — 1, . . . , L— 1, and that of the ML decoder at round L, using the 
union Bhattacharyya bound [40] on a random coded modulation scheme over Q concatenated 
with linear dispersion space-time modulation. For finite T, we obtain similar conditions to those 
in [26]. Finally, as T — > oo, we show that the SNR exponent of random codes is given by the 
Singleton bound for all values of R\ where (|26l) is continuous. ■ 
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Theorem [2] states that optimal diversity gain of N t N r LB and N t N r B for short- and long-term 
models, respectively, can also be achieved by discrete signal sets coupled with linear dispersion 
space-time modulators with constellation Xb (D = 1), namely, space-time modulators that spread 
the symbols of Q over the B fading blocks at each ARQ round. Under this scenario^, full diversity 
is maintained for all rates < R\ < QN t . However, as anticipated in Section Hl-Bl there is one 
drawback of practical concern, namely, complexity. In order to achieve full diversity, the linear 
dispersion space-time modulator needs to spread the symbols of Q over the B blocks, which 
implies that the size of the constellation of each ARQ round is \X B \ = QN t BT. We may, 
however, choose a modulator that spreads symbols over M blocks where M < B in order to 
reduce the complexity of the ML decoder. In this case, there is a tradeoff between the parameters 
of (|26|) . This can be seen as a manifestation of the discrete nature of the input constellation, 
which limits the performance of the outage probability at high SNR. Theorem |2] generalizes the 
result of [22] for the quasi-static MIMO channel to the ARQ block-fading case. 

The upper bound (l26l) is also applicable to any systems using block codes over LB independent 
block-fading periods. The significance of the ARQ framework is that it provides a way of 
achieving the optimal SNR exponent attained by a block code with LB coded blocks, without 
always having to transmit all LB code blocks. Indeed, following [18], observe that 

Pr(^)^Pr(^,...,^ c ) 

< Pr(^ c ) 

= Pr(V>/(y/,H/) = 0) 

< P out (p, £, R x ) + e 

± p -d* D {Ri) (27) 



3 Within our framework, it would also be possible to modulate over 1 < M < BL periods in the short-term case, namely, 
spreading the modulation symbols also across ARQ rounds and l !26t would remain valid. In particular, letting M = LB, we 
could achieve full diversity over the full range of Ri, namely, < Ri < QN t L, which is the same exponent of the Gaussian 
input. However, generalizing our model to this case, would compress key concepts such as information accumulation in a single 
formula, rather than the more natural sum expression in H5\ . In particular, one could define the equivalent channel model up 
to round £ as 

H, = diag (Hi, .... Hi, 0, .... 0) € C LBN " xLBNt . 
where is the zero matrix, and obtain the result. 
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On substitution of (|27l) into (|13l) . we find 

^ L ) ^ ; i ^ (J?l) = fli, (28) 

which shows that the transmit throughput is asymptotically equal to R\ (since R\ > rj(Ri,L)), 
the rate of a single ARQ round. In other words, provided the SNR is sufficiently high, ARQ 
systems which send on average B coded blocks can achieve the same diversity gain as that 
achieved by a block code system which sends LB coded blocks every time. This is because in 
the high SNR regime, most frames can be decoded correctly with high probability based only on 
the first transmitted code block. ARQ retransmissions are used to correct the rare errors which 
occur almost exclusively whenever the channel is in outage. While the throughput r)(Ri,L) is 
a function of L at mid to low SNR, it converges towards R 1 independent of L at sufficiently 
high SNR. Since the optimal diversity gain is an increasing function of L, this behavior can be 
exploited to increase reliability without suffering code rate losses. However, as noted in [18], this 
behavior is exhibited only by decoders capable of near perfect error detection (PED). Therefore, 
the performance of practical error detection schemes can be expected to significantly influence 
the throughput of ARQ systems. 

Since equation (T28l) relates the asymptotic throughput with the coding parameter Ri, the 
optimal SNR exponent given in (T26l) gives the optimal throughput-diversity -delay tradeoff of 
MIMO ARQ block-fading channels^. Examining the optimal throughput-diversity-delay tradeoff 
(l26l) in more detail, we first note that 

Rq Rq _ 
N t LQ ~ QN t ~ T 

is the code rate of a binary code. i.e. < r < 1, as if the coded modulation scheme Cq was 
obtained itself as the concatenation of a binary code of rate r and length N t LQBT. Expression 
(T26l) implies that the higher we set the target rate R\ (equivalently, R ), the lower the achievable 
diversity order. In particular, uncoded sequences (i.e. R\ = QN t L) such as the full diversity 
modulations [37,41], achieve optimal diversity gain of MN t N r , while any code with non-zero 
Ri < QN t L will achieve optimal diversity less than or equal to MN t N r LB or MN t N r B in the 

4 We stress the fact that the coded modulation schemes considered in this paper have a fixed rate, and therefore zero multiplexing 
gain as defined in [17, 18]. However, it is not difficult to show that allowing Q — £logp would imply the achievability of the 
diversity-multiplexing-delay tradeoff of [18]. 
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short- and long-term static models, respectively. This is an intuitively satisfying result as LB 
and B are precisely the number of independent fading periods in the short- and long-term static 
models, respectively, each with inherent diversity MN t N r . 

Figures |2j [3l 0] and [5] are graphs of the tradeoff function (1261) with varying Q, B, L and M 
plotted against the rate of a single ARQ round Ri . We show the tradeoff function ((26)) for both 
short- and long-term static fading models, respectively. 

First we examine the effect of the constellation size Q on the optimal diversity tradeoff 
function. Figure [2] shows the tradeoff curve for three different values of Q. We can see from the 
plot that the tradeoff curves for higher Q are strictly better than lower Q in terms of achievable 
diversity gain. This implies that a high order modulation scheme always outperform lower order 
modulation schemes in the limit of high SNR in terms of error rate performance, for any code 
rate. Alternatively, a system with high Q can choose to operate at higher code rates than a low 
Q system and still maintain the same diversity gain. 

Figure [3] shows the diversity tradeoff curve for different values of B. Similar to the previous 
tradeoff curve with constellation size Q, we observe that systems with high values of B are 
strictly better than systems with low B (in terms of diversity gain). In addition, we notice that B 
corresponds to the number of "steps" in the tradeoff function of (l26l) . Systems with low values 
of B maintain the same diversity gain over wider intervals of rates than systems with high B. 
Relatively, the penalty for using codes with high spectral efficiency is much higher for systems 
with large B (although these systems will still achieve higher diversity gains than systems with 
low B). 

Figure 0] illustrates the effect of the maximum number of allowed ARQ rounds L on the 
diversity of the system. It is clear from the plot that in the short-term static case the effect of 
L is to simply shift tradeoff curves upwards. This is intuitively satisfying, since each additional 
ARQ round represents incremental redundancy, which can be considered as a form of advanced 
repetition coding. Each additional ARQ round contains B additional independent fading blocks 
and hence the diversity gain with L ARQ rounds is simply the diversity gain with L — l rounds 
plus B. On the other hand, in the case of long-term static fading, since the different ARQ rounds 
use the same channel realization, larger L implies a broader range of Ri for which maximum 
diversity can be achieved. 

Figure \5\ shows the impact of M on the tradeoff curve. As anticipated in Section III-BL we 
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observe that the larger M, the larger the optimal SNR exponent. As M increases, larger diversity 
is maintained over a larger range of R\. A careful look to (|26l) reveals that for M > 1, each 
ARQ round behaves as a MEMO block-fading channel with fj = D blocks, each with inherent 
diversity MN t N r , reducing the number of steps of the tradeoff curve. Unfortunately, however, 
increasing M implies an exponential (in M) increase in the overall decoding complexity. 

Remark 1: In [29,42], the authors examined the performance of codes over MIMO block- 
fading channels without ARQ. Using the notation in this paper, the diversity gain based on the 
worst pairwise error rate performance was shown to be upper-bounded by 



The bound in (1291 is based on the fact that the rank of a the codeword difference matrix of a 
given pairwise error event cannot be larger than the minimum number of non-zero rows. The 
application of the Singleton bound [25] to the minimum number of non-zero rows (interpreted as 
the Hamming distance of the code) leads the result shown in (1291 [29,42-44]. Since the bound 
<|29l was derived for the non-ARQ case, we will compare it with our results by letting L = 1 
in (l26l) . An important assumption made in the derivation of (l29l) is that a signal constellation 
of cardinality 2 Q is used for signaling at each transmit antenna. Under this assumption, the 
Singleton bound and the rank criterion give rise to the PEP diversity bound (l29l) . In our case, 
we do not restrict the signals out of each transmit antenna to belong to a constellation of size 
2*3, but rather, allow for more freedom in the system by linearly modulating (combining) MN t T 
-ary symbols to be transmitted over MT channel uses. Figure [6] compares the Singleton bound 
(|29|) with our main result (l26l) . As we see, even in the case of M = 1 our bound yields a larger 
exponent. This effect was also observed in [22] for the quasi-static MIMO channel. 



Having established the main effects of each parameter in (|26|) . we now consider the practical 
coding aspects of Theorem 2. The diversity tradeoff function (|26l) can be viewed as a modified 
version of the Singleton bound [25] with the diversity gain corresponding to the Hamming 
distance of our code C, viewed as a code of length = LD constructed over an alphabet 
of size 2 QMNtT . This is a useful interpretation and naturally leads us to investigate the role of 
Singleton-bound-achieving MDS codes. The role of MDS codes as block codes in block-fading 
channel has been examined extensively in [23,24,26,45]. 




(29) 



VI. Maximum Distance Separable Space-Time Codes 
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In this section, we illustrate that the optimal SNR exponent shown in (|26l) can be achieved with 
practical MDS coding schemes. The block diagram of the concatenated MIMO ARQ transmitter 
structure considered in the numerical examples is shown in Figure [7J A codeword of the MDS 
outer encoder is partitioned into LB blocks. Each such block is then passed through a pseudo- 
random interleaver, subsequently mapped onto a block of complex symbols according to the 
signal constellation, and passed through a linear dispersive modulator. In the ARQ transmitter, 
B blocks of T channel uses are transmitted in each ARQ round. For simplicity, we make use of 
the MDS convolutional codes presented in [23] to illustrate the practical meaning and importance 
of the diversity tradeoff curvej. The ARQ decoder defined in Section III-CI is impractical due to 
the complexity of the typical set decoder. Instead we develop a bounded-distance ARQ decoder 
and a sub-optimal iterative a posteriori probability (APP) based ARQ decoder, respectively, 
approximating the behavior of the typical set decoder. 

For the numerical examples, we consider two systems. The first system has a maximum 
number of ARQ rounds of L = 2, B = 1, and is using the 4-state [5, 7] 8 outer convolutional 
code, while the second system has a maximum number of ARQ rounds of L = 4, B = 1, and 
is using the 4-state [5, 5, 7, 7] 8 outer convolutional code. The rate of the first ARQ round, R\, is 
the same for both systems. The two systems are investigated for both single-input, single-output 
(SISO) and MIMO block-fading channels, subject to short-term static fading and long-term static 
fading, respectively. 

We first consider the use of a bounded-distance ARQ decoder. Define the set of messages 
Vf C X, where the corresponding received codeword hypotheses H^X^(m), m E A4 are within 
a bounded distance from the received matrix Ye, 



Vi=lmeM : 



<£BTNJl + 8)}, (30) 

F 



where 5 > 0. For 1 < I < L — 1, the output of the bounded-distance ARQ decoder is then given 

5 The main goal of these examples is not to approach the outage probability of the channel, but rather to illustrate the meaning 
and significance of the results presented in the previous section. If one wants to approach the outage probability, more powerful 
codes should be employed. For details on outage approaching code ensembles for SISO and MIMO channels the reader is 
referred to [26,46-48]. 
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otherwise 

Denoting the true message rh, the undetected error probability is bounded as 




(31) 



Pr(^,^)=Pr |J (V, = {m» 



(32) 



< Pr (\W e \ 2 F > £BTN r (l + 5)) 

< (1 + 6) eBTNr exp(-£BTN r 5) 



(33) 



(34) 



where (a) follows from bounding the chi-squared distribution of | with the Chernoff bound. 
Finally, letting 5 = /31ogp for (3 > 0, we have 



This result implies that arbitrarily low undetected error probability can be achieved by the new 
decoder, at the cost of additional delay. In particular, (3 should be chosen such that BTN r j3 > 
d*(Ri) in order to achieve the optimal ML exponent d*(Ri). 

Figure [8] illustrates the performance of the two ARQ systems in the short-term SISO static 
channel. We choose the pseudo-random interleaver to be the trivial identity interleaver, i.e. no 
interleaving is applied between the outer encoder and the inner modulator. The mapper over Q 
is set to be BPSK, the space-time modulation rotation matrix R = I, and T = 100 channel uses. 
We apply the list Viterbi decoder proposed in [49] to implement the ARQ decoder outlined in 
(T30l) and (f3TT) . In particular, we choose = d Jj^ t0 minimize the number of retransmissions. 

Considering the L = 2 system, the top three curves in Figure [8] show the corresponding 
outage probability, FER with list decoding and FER with PED. The FER curves are parallel to 
the outage curve at high SNR, which show that the convolutional MDS codes indeed achieve 
the optimal diversity gain. The L = 4 system corresponds to the bottom three curves of Fig. [H 
where again we see that the optimal diversity gain is achieved by the MDS convolutional code. 

Comparing the two ARQ systems, it is clear that significant performance gains can be obtained 
at the expense of higher delays. At FER of 10 -2 , the gain of the L = 4 system over the L = 2 
system is already 5 dB. The performance gap increases even more dramatically at higher SNR. 



BTN r l3 



(35) 
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Figure [9] shows the average number of ARQ rounds of the two ARQ systems considered 
above. For each system, we plot the average number of ARQ rounds with PED, with the list 
decoder and the lower bound given by (fT2)) . respectively. It is clear from the plot that at medium 
to low SNR, significant loss in throughput is incurred by codes that do not approach the outage 
probability limit, like convolutional code. Even more loss in throughput is observed when list 
decoding is used as the error detection mechanism. 

Finally, note that the average ARQ round curves converge towards one at high SNR. This 
agrees with (1281) and shows that regardless of the maximum number of allowed ARQ rounds L, 
no spectral efficiency penalties are incurred at sufficiently high SNR. In the limit of high SNR, 
the transmit throughput r](Ri, L) = R 1 . 

Figure [10] and Figure [TTJ correspond to the error rate and average latency of the same two ARQ 
systems, under long-term static fading. As predicted by the theoretical results of the previous 
section, under long-term static fading both schemes have the same SNR exponent. As a matter 
of fact, despite a 1 dB difference in outage probability, both schemes show virtually the same 
error probability. As already mentioned in the previous section, in the long-term static case, the 
ARQ gain translates in a larger range of R\ supported with optimal SNR exponent. 

We now consider 2x2 MIMO systems with L = 2 and L = 4 using the 4-state [5, 7] 8 and 
[5, 5, 7, 7] 8 convolutional codes, concatenated with the optimal 2x2 linear dispersive modulator 
suggested in [29]. In this example, the channel coherence time is T = 32 channel uses and 
the mapper over Q is set to 4QAM. In this case, the bounded-distance ARQ decoder in (13TI) 
also becomes impractical, and we therefore resort to sub-optimal iterative error detection and 
decoding schemes. As a benchmark, we consider an iterative scheme based on the full-complexity 
APP detector, recursively exchanging code symbol extrinsics with an outer APP decoder, thus 
generating estimates of the information sequence. Applying the max-log APP detector in place 
of the full-complexity APP detector provides a low-complexity alternative. For the examples 
considered here, the full-complexity iterative decoder is roughly twice as complex as the max- 
log APP alternative. For the full-complexity iterative decoder, we only consider PED as the 
target benchmark, while for the max-log APP based iterative decoder we consider PED, as well 
as a non-ideal error detection scheme. At each ARQ round, we run the accumulated received 
signal through six iterations of the respective iterative detection and decoding algorithms before 
examining the decoder output. 
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In the non-ideal error-detection case, errors are detected by examining the soft output of the 
decoder at each ARQ round. Specifically, we use the minimum bit-reliability criterion [50], 
checking at the end of each ARQ decoding round whether the minimum bit-wise log-likelihood 
ratio (LLR) of the information sequence exceeds a threshold, i.e., 

min{|Ly}>0, (36) 

where L e d { denotes the ith element of the information LLR sequence at the £th ARQ round and 
K denotes the length of the LLR vector. If (l36l) holds, decoding is considered successful, and 
the information sequence corresponding to the LLR vector is delivered to the sink. The choice 
of 9 affects both the average latency as well as the error rate of the system. In general, choosing 
a high 9 encourages the receiver to request additional retransmissions, which in turn reduces the 
error rate. However, if 9 is set too high, the system behaves as a block coded system and the 
spectral efficiency advantage of ARQ systems is not realized. Further, it is necessary to increase 
9 as a function of SNR in order to achieve error rate performance comparable to that of perfect 
error detection. To this end, we adjust the threshold as 

9 = max{l, f3 log p}, (37) 

where we have lower bounded 9 in order to encourage retransmissions at low SNR. This choice 
of 9 was found to perform well when the growth parameter (3 is carefully selected. In the 
examples shown here, (3 is determined experimentally. 

Figure [12] compares the error rate performance of the L = 2 system and L = 4 system 
under the short-term fading dynamics. For each system, we plot four curves, corresponding to 
the lower outage probability bound, obtained by using (l24l) . the PED performance for the two 
iterative decoders, as well as the minimum bit-reliability criterion (MinLLR) performance for the 
max-log APP based iterative decoder. In this case we have (3 = 16 and (3 = 32 for the MinLLR 
scheme when L = 2 and L — 4, respectively. We notice that additional retransmissions lead to 
an appreciable decrease in error rates, and, equally important, the MinLLR criterion performs 
virtually as good as perfect error detection. Also, we observe no appreciable loss in performance 
of the max-log APP based iterative decoder as compared to the full-complexity case, confirming 
the use of the max-log APP approximation is well justified. 

Figure [13] compares the average latency (measured in number of ARQ rounds) of the two 
ARQ systems under the short-term fading scenario. Again, we plot four curves per system, 
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corresponding to the lower bound of expected latency, using (1241) and (fl2l) . as well as the PED 
and MinLLR performances. In this case, we observe that the cost of using the MinLLR criterion 
is mainly an increase in latency, caused by requesting superfluous retransmissions, and again 
there is no appreciable loss in performance by applying the max-log APP approximation. 

Figure [14] and Figure \\5\ correspond to the error rate and average latency of the same two 
ARQ systems, under long-term static fading. In this case we have (5 — 12 and (3 = 24 for 
the MinLLR scheme when L = 2 and L = 4, respectively. Once again, as predicted by the 
theoretical results of the previous section, the error rate curves have the same exponent and, 
moreover, have very similar gains. Similarly, the advantage of ARQ in this case is that larger 
throughput can be supported with optimal SNR exponent. 

VII. Conclusions 

The focus of this paper is to derive the optimal tradeoff between throughput, diversity gain, 
and delay for the block-fading MIMO ARQ channel. We prove that for the block-fading MIMO 
ARQ channel with input constellation satisfying a short-term power constraint, the optimal SNR 
exponent is given by N t N r LB for short-term static fading and N t N r B for long-term static fading, 
which is achieved by Gaussian codes of any positive rate. 

When the input signal constellations are constrained to be discrete, this is no longer the case. 
Due to the discrete nature of these signal sets, a tradeoff between rate, diversity and delay arises. 
As our main result, we prove that for the block-fading MEMO ARQ channel with discrete input 
signal constellation of cardinality 2 QNt satisfying a short-term power constraint, the optimal SNR 
exponent is given by a modified Singleton bound, relating all the system parameters. In particular, 
we show that the tradeoff highlights the roles of the complex-plane signal constellation through 
Q, the rate of the first ARQ round R\, the maximum number of ARQ rounds L, and the number 
of fading blocks per ARQ round B. Furthermore, the optimal tradeoff expression includes the 
effect of the space-time spreading dimension M of the linear dispersion modulator, providing 
also a reference of decoding complexity. 

Finally, we present numerical results demonstrating the practical significance of the theoreti- 
cal analysis, showing that practical MDS codes achieve the optimal throughput-diversity-delay 
tradeoff. 
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Appendix 

In this Appendix, we show the details of the proof of Theorem [21 In particular, we detail the 
proof for the short-term static model. The proof corresponding to the long-term static model 
follows exactly the same steps, and it is thus omitted. 

Proof of Theorem [2} Converse 

To prove Theorem [2l we first establish the converse and show that the diversity gain is upper- 
bounded by (T26l) . We assume N t > N r throughout the analysis with no loss in generality^. 

We start following the arguments in [18, Appendix I] and conclude that by Fano's inequality 
we can obtain a lower bound to the error probability of the ARQ decoder at any ARQ round I 
by using an ML decoder that operates over the L ARQ rounds. Therefore, 



Pe (P) > E 



1 I(P\G L ) l__ 

RqL RqLBT 



(38) 



where \x\ + = max{0,x}. Hence, for sufficiently large T, we have that [17, 18] 

Peip^Pon^^Rt). (39) 

Therefore, it follows that we can upper-bound the SNR exponent of the ARQ system by 
considering the outage probability up to ARQ round L. 

Now, we study in more detail the properties of P ou t(p, L, R{) when discrete signal constella- 
tions are used. In particular, we recall that (|24|) states that 

1 D i 1 M ( \ 1 

I{p\Gi) < TJ E min 1 m S log2 det 1 + ^ G ^m+mG[ {d _ 1)M+m (40) 

d=l I m=l \ t / J 

and therefore, 

Pout(P, L, i?i) 

1 D { 1 M ( 

> Pr I > — V min I QN t , — V log, det ( I + -&G a d-i)M +m Gl 



/ l d ( M 

\l=\ d=l L m=l 

(L D ( M N r , 

EE min w.iEE lo M ] 
l=\ d=l l m=l i=l ^ 



-^r^£,(d-l)M+m,i ) } < DR\ 




6 If Nt < N r , it suffices to replace det ^1 + G^bG] by det ^1 + G\ b Gi,b^j in the computation of the input-output mutual 
information with Gaussian inputs and all the arguments still follow. 
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where \e t (d-i)M+m,i < • • • < ^t,(d-i)M+ m ,N r are the ordered A^. eigenvalues of the iV r x 7Y r matrix 
Ge^d-i)M+ m G\ u-x)M+m corresponding to ARQ round £ and fading block (d — 1)M + m. 

We now characterize the behavior of the outage probability at high SNR. Following [17] we 
define the SNR normalized eigenvalues as 



',b i 



logp 

The joint probability distribution of a^j = (a^,&,i, ■ ■ ■ , ott,b,M), can be described using a result 
in [17, Lemma 3] 

N r / N r 



i<j \ i=l 



(44) 



where i^jv t ,/v r is a normalizing constant. Then it follows that 

P out (p,L,i?i) 



>Pr ^^min QJV t , — ]T ]T log 2 I 1 + |-A,, ((J _ 1)M+m , ) [ < (45) 

\£=1 d=l L "1=1 f=l \ * 'J / 

= Pr EE m H ^ m EE lo g^ |1 " a ^ I)M+ra, * l+ < Di M (46) 

\J=1 <Z=1 L m=l i=l J / 



L D C , M TV, 



= Pr(^^nM|QiV t) ^^^|l-a^ 1)JtfM 4|<M 1 ]. (47) 

\£=1 d=l I m=l i=l J J 

If we now define 

aw = «( d -i)M+i, • ■ ■ , «^m)' g R MJVr (48) 

= (a!€,(d-l)M+l,l, • • • j a i,{d-l)M+l,N r , ■ ■ ■ i tt<,dM,l, • • • , Otl,dM,N r ) (49) 

equation (|47T) becomes 

f L d \ 

where a >z b denotes componentwise inequality, i.e., Oj > fej, Vi = 1, . . . , n for some a,b6K n 
and 1 is the all-one vector, since 

f i ogp J£^ 1 (o when a t , d hi 

mm < QN t , — ^— 2^ 2^ I 1 ~ a ^{d-i)M+m, L \ + > = < (51) 
L m=i i=i J QAt otherwise. 
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This means that asymptotically for large SNR, when all the components of ct k}d are larger or 
equal than one (deep fades) the mutual information tends to 0, and to QN t otherwise. Following 
similar steps as in [17] we can write that 



L D M N r 



P on t(p, L, Rx)> / exp - log p y~] V] 5^(2i -l + N t - N r )a e ^ 1)M+m ,i da 



where the large SNR outage event is given by 



1 d=l m=l i=l 



(52) 



6 L = { a e M. LDMNr : E E (! - ^»<^1 (53) 

i=i d=i ^ 1 



a e 



L D 

■>LDMN r 



gp {3wt l } > D ( £ -*L)} 
and a = [a[ 1 , . . . , a L D )' e ^DMN r _ Applying Varadhan's lemma [51] we have that 

( L D M N r \ 

dURi) < _ mf E E E E( 2? - 1 + N * - N r )a Ud - 1)M+ m,i ■ (55) 



OteO L ( 



=1 d=l m=l i=l 

The infimum (1551 ) is solved by considering two cases. If Ri > LQN t , then the infimum is satisfied 
by &e,d = for all I and d, hence the diversity gain is zero. Alternatively, if Ri < LQN t , then 
among all possible vectors a^d, for I = 1, . . . , L and d — 1, . . . , D, we need to have k vectors 
equal to the all-ones vector (a. e d = 1), for some k e Z in order to satisfy the infimum. The 
condition to be met is written 

k > D ( L — ] , (56) 



(57) 



\ QN t/ 

which implies that in order to achieve the infimum k should be 

Ri 



1 + 



D L 



QN t . 

Since ^2 m=1 Yli=i 2i — 1 + N t — N r = MN t N r , we upper-bound the optimal SNR exponent as 



d* D (Ri) < MN t N r [ 1 + 
which proves the desired converse result. 



M V LQN t 
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Proof of Theorem [2j Achievability 

To prove the achievability of the upper-bound on the SNR exponent in (1581) . we examine 
the average frame error rate obtained using random codes and the ARQ decoder described in 
Section Hl-C I This decoder behaves like a typical set decoder for ARQ rounds I — 1, . . . , L— 1, 
and as an ML decoder at round L [38]. Since the channel matrix encompasses the channel 
realizations of all ARQ rounds, with a slight abuse of notation we can express the error probability 
conditioned on the fading realization as 

/ \ 



L-1 



F e (p|H L ) = ^Pr 



2Vi, |J <MY,,H,) 



m 



\ 



J 



+ Pr [Vl..!, |J ^(Y l ,H l ) 



m 



(59) 



where all parameters are defined in Section Hill As shown in [18,38, Appendix I], V5 > and 
sufficiently large T, there exists a code for which the error probability corresponding to the first 
L — 1 rounds can be bounded as 

/ \ 



L-1 

5> 



i=i 



2Vi, U <M Y ^ H 



e =m 



rh^m 



< 5. 



(60) 



Therefore, 



where 



P e (p|H L )<(L-l)5 + P e ml (p|H L 



P e ml (p|H L ) ^ Pr (v L _ u |J ^l(Y l ,H l ) = m J 



(61) 



(62) 



is the error probability of an ML decoding error at the Lth ARQ round. We now characterize the 
behavior of P e ml (p|H^) for a particular code constructiorO, namely, random codes constructed 



7 We could simply conclude the proof by following the same arguments of the proof in [18, Appendix I], namely, using 

p^ l (p\u L )<5+i{u L eO L } 

to argue that P e {p)<p dD ^ Rl ^ (see [18, Appendix I] for details). However, the specific analysis of the ML decoding error 
probability for round L using random codes encompasses the standard quasistatic and block-fading MIMO channels with no 
ARQ as special cases, and therefore is of broader interest. 
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over Q, concatenated with random linear dispersion space-time modulators described in Section 

Following the steps of [22] we consider that the 2 RoLBT codewords of Cq are generated with 
the uniform probability distribution over Q, namely, Vc Q G Cq, 

1 



LBTNj 



p{Cc 



<2| 2Q LBTNt ' 



(63) 



Each codeword Cq G Cq is partitioned into LD vectors, denoted Cg ^ d G qmtni^ w h ere £ 
1, . . . , L and d = 1, . . . , D, such that c Q = [cg^i, . . . , Cq^^]'. Now let 



ft = {R G 



5 MTNt x MTNt 



: RR' = R'R = 1} 



(64) 



denote the set of orthogonal matrices of dimension MTN t x MTNt. As outlined in Section 
III-BL the modulated signals are given by 



Rc 



QAd- 



Then, if we define 



— matMATtxT(x^) 



(65) 



(66) 



where the operator A = mat nxm (a) formats vector a G C nm into an n x m matrix, we have 
that the portion of codeword transmitted over ARQ round t can be written as 



A^i) ■ ■ ■ i-^e 



D 



Then we have that the conditional pairwise error probability is given by 
P (X(n) -> X(fc) \H L = G L )=Q 



< 



exp 



P 
9 

~4N t 



G i (X(n)-X(A;)) 



G L (X(n) - X(fc)) 



It follows from the structure of Gl that 



L D , 

e=i d=i ^ * 



L D 

nri ex p 



M 



||Gf i ( £ i-i)M+m(X^(d_i)Af +m (n) — X^( t j-i)M+m(^)) 



|2 



(67) 

(68) 
(69) 

(70) 
(71) 



If the elements of 1Z are drawn with the uniform probability distribution, it follows from [52, 
Theorem 1] that R has full diversity with probability one, namely, the matrices X^( rf __i) A / +m (n) — 
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~^-e,{d~i)M+m(k) have full ranlc|. We now apply the singular value decomposition (SVD) [36] to 
both channel and difference matrices 



and 



^e,(d-i)M+ m (n) - y^z,(d~i)M+m{k) - AD| (d _ 1)A/+m B t 



(72) 



(73) 



and get that 

P fx(n) -> X(fc) | H L = G L 



L D 

<nri ex p 

1=1 d=l 
L D 

=nri ex p 

l=\ d=l 
L D 



M 



AN, 



i,(d-l)M+m\X-l,{d-l)M+m 



(n) - X £ 

l)M+m 

(*)) 



|2 



P 



m=l 
M 



nn 



exp 



=1 d=l 



P 



E ll UA l 



V1"A"D 2 Rt 



m=l 
M 



1 1 

2 pf)2 

^(d-lJM+rn ^ t,{d-l)M+m 



m=l 



(74) 
(75) 
(76) 



i i 

where P = V+A is unitary. The diagonal matrices A|, d _^ M+m and D| ( - d l ^ )M+m are composed 
of the singular values of the channel matrix G^(d-i)M+m and codeword difference matrix 
^e,(d-i)M+ m (n)-^e,(d-i)M+ m (k), respectively. As mentioned earlier, the matrices X^ i((i _ 1)M+m (ra) 



X 



i,(d-l)M+m 



(k) have full rank with probability one, which implies that the the MN r singular 



values in D| (d _ 1)Jvf+m 
we now define 



are all non-zero for m = 1, . . . , M, d = 1, . . . , D and t = 1, . . . , L. If 



e,(d-i)M+m 



PD 



t,(d-l)M+m 



Pt 



and 



ll,{d-l)M+m - dia § (Te,(d-l)M+m) 

— (le,(d-l)M+m,li ■ ■ • > lt,{d-l)M+m,Nr)i 



(77) 

(78) 
(79) 



8 As it will be clear in the following, random rotations are not essential in the proof. It is sufficient to rely on the existence 
of a particular R with full diversity [29, 32-35, 37, 41]. 
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we can rewrite (1761) as 



L D 



M N r 



P (X(n) - X(fc) | H L = G L < H H ex P " E E 



=1 d=l 



m=l j=l 



7l,(d-l)M+m,i l-a £ , (d _ 1)M+m , I 

m 1 



Averaging (l80l) over the code ensemble, namely X(n),X(/c) and R, we get that 



L D 



F(X(„).X W |a) = niI^ 



l=\ d=l 



2QMTN t 



+ 



1 



2QMTN t 



E 



E 



R 



c Q/,d(")^CQ,«,d(fc) 



ex p - E E 



m=l i=l 



7l,(d-l)M+m,i l-a f , (d _ 1)M+m , I 



If we now sum over the 2 RoBLT codewords, we have the union bound 



L D 



=1 d=l 



2QMTN t 



+ 



2QMTN t 



E 



E 



R 



c Q,e,d( n )^ c Q,e,d( k ) 



exp 



M JVr 

EE 

m=l i=l 



7l,(d-l)M+m,i l-« £i(d _ 1)M+m ,, 
4iV, ' 



+ 



exp | -LDMTQN t log(2) 
1 



i? 

QA^ ~ Tdmtqn, 



L D 



EE lo & i 



=1 d=l 



2QMTN t 



E 



E 



R 



c Q,e,d( n )^ c Q,e,d( k ) 



exp 



M N r 

EE 

m,=l i=l 



7<,(d-l)Af +m,i l-a Ud _ 1)M+m:i 

AN 1 



= exp (-LDMTQN t hg(2)E(p, a)) 
where we have defined the union bound exponent as 

L D 



S(p,a) = 1 - 



r T(l N. 



QiV t LDMTQN t j^j^ 



+ 



1 



2QMTN t 



E 



E 



R 



c s,«,<i( n )/ c aA<iW 



ex p - E E 



m=l i=l 



7€,(d-l)M+m,i l-a,. (d _ 1)M+ma 



(80) 



(81) 



(82) 
(83) 



(84) 
(85) 

(86) 



(87) 



Following similar arguments to those in [22,26], we use the dominated convergence theorem 
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[53] to obtain that 

M N r 

lt,(d-l)M+m,i i 



lim E R 

p^OO 



exp ~ Yl Yl m,i p l_ ^' (d - i)M+m ' 1 



AN t 

m=l i=l 



E 



R 



p^co 



M N r 

Jim exp ( YY''"' 1 ix""'' l y ' UV 



(88) 

m=l i=l / . 

= 1 - i{a M 5- 1}, (89) 

since ^(d-^M+mj > with probability one. For e > and large SNR, the union bound exponent 
E(p, a) can be lower-bounded by 

^ * £=1 d=l 

Let now 

£ e = {a G R LDM ^ : £ e (p, a) < 0} (91) 
= |aGR^ : EEl{^^l-e}>LD(l-^|. (92) 
Then we can bound the overall error probability as 

P e (p)< I min{l, exp (-LDMTQN t \og(2)E e (p, a))} p(a)da. (93) 

In a similar way to what it is done in [26], we consider codes with block length T(p) such that 

r ± lim pBl. (94) 

p^oo log p 

That is, we consider sufficiently long codewords large SNR such that the error probability is 
never dominated by the event when two codewords coincide. Thus, we can write that, 

L D M N r 



P e (p)< / exp - logpV V V V(2i - 1 + JV t - iV r )a £i((i _ 1)M+mii da 

J(X£tenR + \ e=ld=lm=li=1 J 

- / L D M N r 

+ / exp -logpVVVV(2z-l + iV t -7V r )a, )l 



,(d-l)M+m,i 

'ae£ ? rw— \ ~ ~ ~i ~ 

+ rLDMQNt log(2)£ e (p, a) j da (95) 
and therefore, the random coding exponent is lower-bounded by 

d (r) > sup min{ d u d 2 } (96) 
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where 

~ L DM N r 



d * = in L MN E E E E( 2 * - 1 + ^ - iV >'.^-i)«+™^ ( 97 ) 



+ ^ 1=1 d=l m=l i=l 



is the exponent for large enough codewords and 

{L D M N r 
E E E E^ " 1 + Nt ~ N r)^(d-l)M+m,i 

+ TLDMQN t \og(2)E t (p,<x)^ (98) 

inf I rLDMQNt log(2) ( 1 - — - 
aeesr&l DUNr { V Q N t. 

L D \ 

+ M(N t N r - rQN t log(2)) ^ ^ l{a M ^ 1 - e} I (99) 

is the exponent that characterizes the finite block length. 

Following similar steps to those in the converse, the SNR exponent of the first component d\ 
can be written 



di > (1 - e)MN t N r 



LB f i R x 



(100) 



M \ LQN t/ 

Following similar arguments as in [26], we see that if < rQN t log(2) < N t N r then the infimum 
(l99l) is given by 

LDMrQN t log(2) (l - . (101) 

Otherwise, if rQN t log(2) > N t N r , then the infimum is 

Rn 



TLDMQN t log(2)(l-^)+M(N t N r (l-e)-rQN t log(2))[ LD [ 1 - ^ | -1 

(102) 

The random coding SNR exponent lower-bound can be tightened by letting e — > 0. By collecting 
the results together, we see that for sufficiently large r, ri 2 coincides with d\. In fact, one observes 
that for T — > oo, the overall error probability is given by the probability of the event S e , since the 
second integral in (|951) vanishes. Hence the diversity lower-bound coincides with the diversity 
upper-bound (|26|) for all rates except at the discontinuities. 
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Fig. 1. MIMO ARQ system model. 
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Fig. 2. Optimal diversity tradeoff curve corresponding to L = 2, B — 4, M = 1 for a 2 x 2 MIMO channel. 
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Fig. 3. Optimal diversity tradeoff curve corresponding to L = 2, Q = 2, M = 1 for a 2 x 2 MIMO channel. 
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Fig. 4. Optimal diversity tradeoff curve corresponding to B = 2, Q = 2, M = 1 for a 2 x 2 MIMO channel. 
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Fig. 5. Optimal diversity tradeoff curve corresponding to B — 4, Q = 2, L = 2 for a 2 x 2 MIMO channel. The curves with 
8 steps correspond to M = 1, those with 4 to M = 2 and those with 2 to M = B — 2. 
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Fig. 6. Comparison of the optimal diversity tradeoff curve for with J29b (dashed steps) for B — 4, Q = 2, L = 1 and M = 1, 4 
in a 2 x 2 MIMO channel. 
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Fig. 7. Block diagram of the concatenated MIMO ARQ architecture. The interleaver corresponding to ARQ round t and fading 
block b is denoted by ty^^. 
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Fig. 8. FER with MDS convolutional code over a short-term static SISO channel corresponding to B = 1, Q = 1 and T = 100. 
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Fig. 9. Average number of ARQ rounds for MDS convolutional codes over a short-term static SISO channel corresponding to 
B = 1, Q = 1 and T = 100. 
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Fig. 10. FER with MDS convolutional code over a long-term static SISO channel corresponding to B = 1, Q = 1 and 
T = 100. 
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Fig. 11. Average number of ARQ rounds for MDS convolutional codes over a long-term static SISO channel corresponding 
to B = 1, Q = 1 and T = 100. 
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Fig. 12. FER with MDS convolutional code over a short-term static 2x2 MIMO channel corresponding to B — 1, Q — 2 and 
T — 32. The thick solid lines are the lower outage probability bounds. For L = 2, diamonds correspond to full-complexity APP 
detection with PED, while squares and crosses correspond to max-log APP detection with PED and MinLLR, respectively. For 
L = 4, pentagrams correspond to full-complexity APP detection with PED, while circles and asterisks correspond to max-log 
APP detection with PED and MinLLR, respectively. 
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Fig. 13. Average number of ARQ rounds for MDS convolutional codes over a short-term static 2x2 MIMO channel 
corresponding to B = 1, Q = 2 and T = 32. The thick solid lines are the lower bounds on expected latency. For L — 2, 
diamonds correspond to full-complexity APP detection with PED, while squares and crosses correspond to max-log APP detection 
with PED and MinLLR, respectively. For L = 4, pentagrams correspond to full-complexity APP detection with PED, while 
circles and asterisks correspond to max-log APP detection with PED and MinLLR, respectively. 
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Fig. 14. FER with MDS convolutional code over a long-term static 2x2 MIMO channel corresponding to B = 1, Q = 2 and 
T = 32. The thick solid lines are the lower outage probability bounds. For L = 2, diamonds correspond to full-complexity APP 
detection with PED, while squares and crosses correspond to max-log APP detection with PED and MinLLR, respectively. For 
L = 4, pentagrams correspond to full-complexity APP detection with PED, while circles and asterisks correspond to max-log 
APP detection with PED and MinLLR, respectively. 
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Fig. 15. Average number of ARQ rounds for MDS convolutional codes over a long-term static 2x2 MIMO channel 
corresponding to B = 1, Q = 2 and T — 32. The thick solid lines are the lower bounds on expected latency. For L = 2, 
diamonds correspond to full-complexity APP detection with PED, while squares and crosses correspond to max-log APP detection 
with PED and MinLLR, respectively. For L — 4, pentagrams correspond to full-complexity APP detection with PED, while 
circles and asterisks correspond to max-log APP detection with PED and MinLLR, respectively. 
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